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SUMMARY 


The  objective  of  this  research  is  to  develop  practical  design  proce- 
dures  that  can  be  used  in  conjunction  with  optimal  digital  control¬ 
lers  and  estimators  for  future  orbiting  large  space  structural 
systems.  In  practice,  observational  data  used  to  verify  the  orienta¬ 
tion  and  shape  of  large  flexible  systems  will,  in  general,  be  col¬ 
lected  on  a  sampled  basis  (discrete-time  data  system).  Random  noise 
is  also  included  with  such  observational  data.  Systems  which  will 
be  designed  to  control  both  the  overall  orientation  as  well  as  the 
shape  of  some  of  the  subsystems  (such  as  an  antenna  mesh  form)  must 
function  in  the  discrete-time  domain.  The  aim  of  this  research  is 
to  develop  methods  that  can  be  used  to  design  linear  quadratic  regu¬ 
lator  (LQR)  or  linear  quadratic  Gaussian  (LQG)  controllers/ 
estimators  for  orbiting  large  flexible  systems.  The  digital  control 
of  an  orbiting  shallow  shell  antenna-type  system,  as  well  as  flexi¬ 
ble  platform- type  systems  are  examined  in  this  report. 
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INTRODUCTION 


The  purpose  of  this  research  Is  to  study  and  develop  the  practical 
methods  of  designing  LQR  and  LQG  digital  controllers  and  estimators 
for  orbiting  large  space  flexible  structural  systems. 

A  number  of  1 nvest i gators E 1 3  have  considered  the  problem  of 
developing  optimal  control  laws  for  large  flexible  orbiting  systems 
under  the  assumption  that  the  state  vector  information  would  be  ob¬ 
served  directly  or  that  the  state  information  could  be  estimated  on 
a  continuous  basis.  In  practice,  however,  the  observational  data, 
in  general,  will  most  probably  be  collected  on  a  sampled  basis  i.e. 
a  discrete-time  data  sys tem . E 7 - 1 0 J  Sometimes  the  random  noise  which 
is  included  in  the  observational  data  cannot  be  neglected.  In  this 
case,  the  observational  data  must  be  treated  as  observational  data 
with  random  measurement  noise.  Therefore,  it  appears  useful  and 
timely  to  study  the  control  problem  of  orbiting  large  space  structur 
al  systems  with  discrete- time  observational  data  and  random  measure¬ 
ment  noi se  .  t 1 1 •  12] 

The  analysis  and  design  of  LQR  digital  controllers  and  LQG  digital 
controllers  and  estimators  have  been  finished  for  a  hypothetical 
thin  free-free  long  beam  in  orbit.  The  concepts  of  the  degree  of 
controllability  (observability)  which  is  a  function  of  the  controlla 
bility  (observability)  matrix  L 1 3  —  1 6 ]  have  been  applied  to  solve  the 
problem  of  actuator  and  sensor  placement.  The  closed-loop  dynamics 
of  the  LQR  and  LQG  control  system  for  a  long,  slender  orbiting  flexi 
ble  beam  have  ben  simulated.  The  design  of  the  digital  controllers 
and  estimators  has  been  certified  by  simulation. 

The  aim  of  this  reported  research,  for  the  period  July  1988  -  July 
1989,  is  to  develop  methods  that  can  be  used  to  design  linear  quad¬ 
ratic  regulator  (LQR)  or  linear  quadratic  Gaussian  (LQG)  control¬ 
lers/estimators  for  orbiting  space  flexible  antenna/ref lector  and 
platform  systems.  The  specific  tasks  reported  are  as  follows: 

(1)  The  analysis  and  design  of  the  LQR  and  LQG  digital  controllers 
and  estimators  for  large  space  flexible  shallow  spherical  shell 
systems  representat i ve  of  orbiting  r efl ector /antenna  systems.  A 
paper  based  on  thi3  task  has  been  accepted  for  presentation  at  the 
40th  International  A stronaut i ca 1  Congress  and  is  the  basis  for 
Chapter  2  of  this  report. 

(2)  The  development  of  three  definitions  of  the  degree  of  controlla 
bility  (observability)  which  are  based  on  the  scalar  measure  of  the 
controllability  (observability)  matrix  for  discrete- time  systems. 
Their  general  properties,  together  with  the  simple  physical  and  geo¬ 
metrical  1 nterpretati ons  for  the  fuel  optimal  control  problem  are 
shown.  The  applications  of  these  definitions  for  the  actuator  place 
ment  problem  of  the  orbiting  shallow  spherical  shell  system  are 
illustrated  for  seven  sample  cases.  A  paper  based  on  this  task  has 
been  presented  at  the  12th  Biennial  Conference  on  Mechanical  Vibra¬ 
tion  and  Noise  and  is  the  nucleus  for  Chapter  3  of  this  report. 


(3)  Modelling  of  large  flexible  orbiting  platform  control  systems 
and  the  analysis  and  design  of  LQR  digital  controllers  for  large 
flexible  orbiting  platform  systems.  This  is  in  progress  and  will  be 
the  basis  of  a  Master's  thesis,  currently  under  preparation.  A 
brief  summary  of  this  task  is  presented  in  Chapter  4  of  this  report. 
It  is  intended  that  the  Master's  thesis  will  also  be  published  as  a 
supplement  to  this  volume. 

(4)  Finally,  Chapter  5  will  summarize  some  concluding  statements 
and  follow“on  plans  for  the  continuation  of  this  general  area  of 
research  along  t  .e  lines  of  our  recent  proposal^'}  to  WRDC. 
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2 


THE  OPTIMAL  LQG  DIGITAL  SHAPE  AND  ORIENTATION  CONTROL  OF  AN 
ORBITING  SHALLOW  SPHERICAL  SHELL  SYSTEM 


2.1  Introduction 


Future  proposed  space  missions  would  involve  large  inherently 
flexible  systems  for  use  in  communi cations ,  radiometry,  and  in 
electronic  orbital  based  mail  systems.  The  use  of  very  large 
shallow  dish  type  structures  to  be  employed  as  receivers/reflectors 
for  these  missions  has  been  suggested.  In  order  to  satisfy  mission 
requirements  control  of  the  shape  as  well  as  the  overall  orientation 
will  often  be  required.  The  proposed  paper  is  devoted  to  a  study  of 
the  LQG  digital  optimal  control  of  the  shape  and  orientation  for  an 
orbiting  shallow  spherical  shell  system. 


The  dynamical  system  equations  and  deterministic  control  of  the 
three-dimensional  shallow  spherical  shell  in  orbit,  in  which  the 
yaw,  roll  and  pitch  attitude  angles  of  the  undeformed  shell  and  six 
flexible  modal  parameters  are  modeled,  have  been  studied.^*  18“19] 
In  this  paper  the  LQG  digital  optimal  control  technique  will  be 
applied  to  the  shape  and  attitude  control  of  shallow  spherical  shell 
systems.  In  order  to  realize  estimation  and  control,  the 
mathematical  models  for  the  actuators  and  sensors  are  developed,  in 
which  up  to  12  jet  actuators  are  used  for  control  of  the  attitude 
and  shape  of  the  shell.  For  observation,  two  sun  sensors,  two  earth 
horizon  sensors  and  six  displacement  sensors  are  assumed  to  be  used 
for  measuring  the  attitude  and  displacement  of  the  deformed  shell. 


As  for  the  suboptimal  digital  LQG  problem^12!  we  are  not  only 
interested  in  the  separate  design  of  both  the  controller  and 
estimator,  but  more  interested  in  determining  the  influence  of  the 
different  combinations  of  the  controller  and  observer  pole  locations 
on  the  estimate  and  control  process  and  the  best  combination  of 
controller  and  observer  pole  locations. 


In  addition,  the  problem  of  actuator  placement  is  discussed  by  means 
of  the  concept  of  the  degree  of  controllability.^20^  The  best 
arrangement  of  the  locations  for  the  actuators  will  be  discussed. 

Finally,  the  simulations  certify  the  analysis  and  design  of  the 
digital  optimal  controllers  and  estimations. 

2.2  Mathematical  Model 


2.2.1  Dynamical  Equations 

The  mathematical  model  of  an  isotropic  shallow  flexible  spherical 
shell  in  orbit  was  developed  in  references  Cl*  18-19],  The  main 
assumption  is  that  the  shell’s  elastic  displacements  were 
principally  in  the  transverse  direction  (parallel  to  the  symmetry 
axis)  and  were  small  as  compared  with  the  other  characteristic 
dimensions  of  the  shell.  The  assumption  of  shallowness  further 
insures  that  the  ratio  of  the  displacement  of  the  shell's  apex  point 
above  its  base  plane  (H)  is  small  as  compared  with  the  radius  in  the 
base  plane,  l,  (Figure  1). 
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The  resulting  linearized  equations  of  motion  for  the  rigid 
rotational  and  generic  elastic  modes  were  developed  as: 

T-ui f-( 1 )J-Cx/Jx°^«c 

**  e  f  Q  )  ? 

4>  +  1»U3<>+(  I-U3)  ?-Cz/J2  «c 

0-3u20-(2t/J^O)  )l  I1(n)en-Cy/Jy0)uo  (D 

n*  1 

ei+(Qi-3)ei+(2I1(i)/Mit)0-3Ii  1  } /Mi  1  +  E  ^Mi  tu*  (1  -  1,  2,  ...  6) 
where 

t  ■  w q t  ,  ej  ■  qj(t)/t  (i  ■  1t  2,  ...  6) 

The  derivative  in  Equation  (1)  is  with  respect  to  x. 

<p,  <p,  9  yaw,  roll  and  pitch  angles,  respectively,  between 

the  undeformed  axes  of  the  shallow  spherical  shell 
and  the  axes  of  the  orbiting  local  vertical/local 
horizontal  system 

q^(t)  modal  amplitude  of  the  ith  generic  mode 

<uc  orbital  angular  rate,  constant  for  assumed 

circular  orbit 

characteristic  length  (the  base  radius) 
the  ith  modal  mass 

principal  moments  of  inertia  of  the  undeformed 
shell 


C  X  •  ^y*  C  z 

*c 

IM  ,  U2  .  W3 


the  components  of  external  torques 

coordinate  of  differential  area  on  the  surface 
above  the  base  plane 


(40)-40))/40).  (40>-40,)/j'0) 


respectively 


4U 

<px 

wn 


transverse  component  of  the  ith  modal  shape 
function 


natural  frequency  of  nth  mode 
5 


un/«Jc 


nn 

The  frequencies,  u ,  of  the  spherical  shell  are  evaluated  usina  the 
following  identities,  as  presented  by  Johnson  and  Rei saner . L2i ] 

wpj  "  ^hTT^  UPJ  +  w-  (2) 

where 

D  -  Eh3/i 2( 1 -D2 )  ,  ending  stiffness  factor 

o>2  -  E/pR2 
00 


v  -  Poisson's  ratio 


p  -  density  of  shell  material 
E  -  modulus  of  elasticity 
h  -  wall  thickness  of  the  shell 
1  -  the  base  radius  of  the  shell 

R  radius  of  curvature  of  the  shell 


The  ppj  are  calculated  from 


14* 


Ip+i’ Upj 


1-v  for  p-0,1 


(3) 


For  the  sample  calculations  in  this  paper,  we  will  consider  only 
such  values  of  ppj  (j  -  1,  2,  3)  for  the  cases  where  p  -  0,  1. 

The  mode  shapes  of  the  transverse  vibrations  of  a  shallow  spherical 
shell  with  a  completely  free  edge  are  given  by:^21,  22] 

4n)  *ApJ  {  cpj  £P  +  Jp  (  PpJ  »  5  )  +DpJ  ^  Pp  j  *  ^  )}  COS  p(B  +  8o) 

pJ 


where  p  the  number  of  nodal  diameters  (meridians) 

J  the  number  of  nodal  circles 

For  the  following  proposal  numerical  values  considered:  H  -  1m;  t  - 
100m:  h  -  0.01m;  v  -  1/3;  R  -  5000m;  p  -  31  .3309  kg/m3  ;  E  -  0.840466 

x  1o9  Newton/m2;  the  frequencies  and  the  shape  function  coefficients 

are  shown  in  Table  1.  (Cpj  -  0;  p  -  0,  1;  J  -  1,  2,  3) 
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TABLE  1.  FREQUENCIES  AND  COEFFICIENTS  OF 
SHAPE  FUNCTION 

"  ^ 

n(p,J)  p  J  u2  (rad. /a)  Apj  Dpj 

PJ 


0 

1 

9.076 

1 

.02778 

2.1979  -0.8381 xIO*1 

0 

2 

38.507 

1 

.02946 

3.1389  0. 3119x10“2 

0 

3 

87.82 

1 

.03692 

3.8468  -0.12787x10*3 

1 

1 

20.52 

1 

.02163 

3.8359  -0. 1901x10*’ 

1 

2 

59.86 

1 

.031 98 

4.9627  0.7045x10*3 

1 

3 

119.00 

1 

.04459 

5.8752  -0. 2845x10** 

2.2.2  Point  Actuator  Model 

The  point  actuators  are  modelled  as  follows.  Let  the  ith  actuator 
be  located  at  (txi>  fcyi.  tZi )  anil  the  force  component  of  the  ith 
actuator  be  (fxj,  f Ui  *  ^zi^ffil  •  The  torques  provided  by  the  a  actu* 
ators  are  as  follows: 

/  C  x  \  U  /  0  *  i.  1  iy  i  \  H  X  i\ 

c  -  Uy )  -  M  *zi  0  4xl  fyl  |  f  i|  (5) 

\Czy  i  ^xi  0  J\?zlJ 

The  corresponding  generic  force  in  the  nth  mode, 

En  -  -4n)(Ift/M)  6 ( x  “xi»  y*Yi»  z-ZiJdm  (6) 

i 

For  the  shallow  spherical  shell  it  is  assumed  that  the  major  elastic 
displacement  occurs  in  a  direction  normal  to  the  base  (y,  z) 

plane,  i.e.,  4>xn^(£,  B)i.  The  control  forces  provided  by  the 

a  actuators  are 

f  =  If i  -  l  (fXi  i  +  fyl  3“  ♦  fzi  k)  |  f i|  (7) 

i  i 

Substituting  (7)  into  Equation  (6),  we  have 

En  -  I  ^n)(^i.  6i)fxi|  f i|  (8) 

1 

2.2.3  Point  Sensors  and  the  Observational  Model 

In  order  to  describe  the  geometric  relationship  between  the  orienta-* 
tion  and  reference  systems,  it  is  necessary  to  define  the  following 
coordinate  reference  systems: 

i$i,  Ji»  f<i*  E}  Earth  inertial  coordinate  system,  its  center  is  the 
center  of  Earth,  “i  i  axis  points  to  the  spring  equinox. 

{iQ*  ^o»  kQ,  o}  Orbital  reference  system,  its  center  is  the  mass 
center  of  the  undeforraed  shell,  the  is  along  the  outward  local 
vertical  from  the  center  of  the  Earth. 

1 1 ,  J,  £,  o}  The  body  reference  system  of  the  undeformed  shell, 
i.e.,  the  principal  axes  system  of  the  undeformed  shell. 

The  relationships  between  the  coordinate  reference  systems  are  as 
follows: 


8 


0)1  Ji 


where  up  is  the  argument  of  perigee,  fe  is  the  true  anomaly  angle,  i 
is  the  orbital  inclination,  and  Q  is  the  right  ascension  of  the 
ascending  node  line. 


J  -  Ab0U,  0,  ♦>  J 


where 


CO34CO30  si  n$cosiJ>  +  co3$3i  n0si  mp 


Abo  ■  **3in*cos0  cos<>cosi(i-sln(tsin0sim|» 


si  n0 


-cos0si  nip 


si  n$si  n4»-cos*si  n0cos<|/ 
cos$si n*+3in$sl n0cos* 


COS  0COS  Ip 


It  is  assumed  that  two  earth  horizon  sensors  are  used  for  measuring 
the  angle,  0e2»  between  the  axis  J  and  t0,  and  the  angle,  @43,  be** 
tween  the  axis  £  and  iQ.  The  observational  geometrical  relation-* 
ships  are  as  follows: 


cos0e2  -5  *  tQ  ■  -  sin$  cos0 
cos0e3  -  £  •  fQ  -  sin0 


(11) 


It  is  assumed  that  two^sun  sensors  are  used  for  measuring  the  angle, 
0S2 ,  between  the  axis  J  and  the  direction  vector  of  the  sun,  $,  and 
the  angle,  ©33,  between  the  axis  £  and  The  observational  geo¬ 

metrical  relationships  are  as  follows: 

cos032  *  Jo*®  ■  “  3in$co30(to«s)  +  (co3$cosi|>-3in$3in03in*)  (J0»s) 

♦  (cos<|>3imp+3in<|>3in<|>sin0co3ip)(£o*3)  (13) 

cos033  -  £*S  -  3lne(fo*§)-cos0sini(i(fo*$)«ico30cosi|>(£o*$ )  (14) 

It  is  assumed  that  Y  displacement  sensors  are  used  to  measure  the 
shell's  transverse  displacement  parallel  to  the  x  axis,  i.e., 

gi  •  CfU ($i ,  8 i ;  t)  •  i 

-  Cf  l  ^J)(5i,  6  i  )qj  ( t )  (15) 

J 


r  .  r  i  J  / 

Cf*  l  *x  (*i.  Si )  e j ( t ) 


(i  -  1  ,  2 . v) 
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where 


U(CiiSi;t)  displacement  vector  of  the  shell  at  (£i.  Si) 

Cf  coefficient  relating  the  displacement  to  the  output  volt* 

age 

l  the  base  radius  of  the  shell 

The  linearized  version  of  Equations  11-14  and  Equation  15  are 
combined  into  the  observation  equation  for  the  shell  system  as 
follows : 


Y 

where 


Hx 


^  CO30e2 

COS!e3  A, 

COS0a2-(jO*3) 
cos0S3-  ^o*§) 
81 
82 


(16) 


-1 

o 


(k0’3)  -(T0*a) 
-(Jo*5) 


(io*S) 


0„3 


1*3 


Cft 


O' 

v-/4xfi 


(1) 


(2) 


(6) 


(♦x  (Cii8i)*x  ^ i • ® i a 

<U  (2)  (6) 

(♦x  U2.S2)*X  U2.S2),...*x  U2,S2>| 


O)  (2)  (6) 

(♦x  ,8  )*x  (C  .8  ),...*x  (C  ,8  )| 

Y  Y  Y  Y  Y  Y 


Jtx6 


0 

W(#*Y)x9 


2.2.4  The  State  Equations  of  the  Shallow  Spherical  Shell  System 

Considering  Equations  5  and  8,  Equation  '1  can  be  written  in  the 
state  vector  form: 

x  -  Ax  +  Bu  (17) 
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where 


x  *  (♦»  ♦*  0.#  €  j  »  ®  2  *  ® 3 •  •••  ®6»  'fr »  ^  •  ® •  ®  1  •  *2  •  •••  ®6 ^ * 

u  •  ( I  f  1 1  •  I  ^  •  *••  I  ^  al  ^ 


The  state  Equation  (17)  can  be  discretized  as  follows 

x  (k*1 )  -  e(AT)x(k)-*T(AT)u(k)  < 1 8 ) 

where 

A  AT 

♦(AT)  -  e  ;  AT  sampling  time 

fAT 

r ( AT )  -  ♦(AT-t)  B  dt 

j  0 

The  discretized  observation  equation  is  as  follows 
Y (k )  -  Hx(k  )  (19) 

2.3*  The  Placement  of  Actuators 

The  concept  of  the  degree  of  controllability  has  been  used  for  deter 
mining  the  placement  of  actuators  on  the  shell.  The  detailed  re** 
suits  of  this  study  are  described  in  reference  [20].  In  this  paper 
we  only  list  part  of  the  main  results  in  Table  2. 

The  degree  of  controllability  is  the  scalar  measure  Of  system 
controllability  and  its  reciprocal  indicates  the  control  effort  of 
the  system.  The  degree  of  controllability,  ui  ,  in  Table  2,  is 
defined  as  follows: 

U1  *  *min^c J 
where 

wc  -  PCPCT 

pc  -  (r  er  e2r  ...  en-1r) 


1 1 


(20) 


TABLE  2.  THE  DEGREE  OF  CONTROLLABILITY  FOR  DIFFERENT  LOCATIONS  OF 
ACTUATORS  ON  THE  SHALLOW  SPHERICAL  SHELL 


Case  No.  Actuators  Locations  of  Actuators  Degree  of  Controllability,  ui 


0.0 

0.13875x10‘12 
0.1 4332x10“ 10 

0. 1 4332x1 0*10 

0.28673x10“10 

0.85391x10“12 

0.45342x10*10 

0.45349x10“10 


The  values  of  the  degree  of  controllability  for  each  case  are  listed 
in  Table  2.  The  reason  why  the  degree  of  controllability  for  Case  1 
is  zero  is  that  all  the  actuators  are  located  along  the  meridional 
nodal  line  of  one  of  the  fundamental  shell  vibrational  modes.  The 
reason  why  the  degree  of  controllability  for  Case  6  is  so  low  is 
that  the  actuators  are  all  located  on  the  same  nodal  circle  (£  - 
const.) . 

It  is  assumed  that  the  thrusters  which  are  located  at  the  shell's 
edge  have  two  Jet  directions:  (i)  tangent  to  the  shell's  edge  and 
(ii)  normal  to  the  shell  surface.  Each  of  the  other  thrusters  has 
only  one  jet  direction  -  normal  to  the  shell's  surface. 

If  the  number  of  actuators  is  limited  to  12,  the  arrangement  of  actu¬ 
ators  as  in  Case  7  or  Case  8  is  suggested. 

2.4  Analysis  and  Design  of  the  LQG  Optimal  Regulators  and  Ob¬ 
servers  for  the  Shallow  Spherical  Shell  System 

The  dynamic  equations  and  observational  equations  of  the  shallow 
spherical  shell  system,  as  a  deterministic  model,  have  been  de¬ 
veloped  in  the  previous  section.  They  are  Equations  18  and  19.  In 
fact,  dynamic  systems  are  driven  not  only  by  our  own  control  input, 
but  also  by  disturbances  which  we  can  neither  control  nor  model 
deterministically.  Sensors  generally  do  not  provide  exact  readings 
of  desired  quantities,  but  Introduce  their  own  system  dynamics  and 
distortions  as  well.  Furthermore,  these  devices  are  also  always 
noise  corrupted.  In  order  to  solve  these  difficulties,  the  theory 
and  method  of  the  LQG  problem  (design  of  the  optimal  stochastic  con¬ 
troller  for  a  problem  described  in  terms  of  linear  system  models, 
quadratic  cost  criteria,  and  Gaussian  noise  models)  are  effective. 

In  this  section  the  LQG  theory  is  applied  to  the  design  of  sub- 
optimal  regulators  and  estimators  for  the  shallow  spherical  shell 
system  . 


2.4.1  Model,  Problem  and  SolutlonC1^] 

Let  a  system  be  adequately  described  by  the  n-dimensional  state 
difference  equations: 

x  ( k  + 1  )  -  $x(k)  +  T(k)  +  w  ( k  )  (21) 

where  the  p-dimensional  u(k)  is  the  control  input  to  be  applied  and 
w(k)  is  the  n-dlmensional  zero-mean  white  Gaussian  discrete- time 
noise  with 

E ( w ( i  )  wT(j  )  }  -  Q  « (22) 

and  assumed  to  be  independent  of  the  initial  state  condition.  The 
Initial  x(0)  is  modeled  as  a  Gaussian  random  vector  with  mean  x(0). 
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and  covariance,  Pq .  Available  from  the  system  are  m-dimensional 
sampled-data  measurements  of  the  form 

Y (k )  -  Hx (k  )  ♦  v (k  )  (23) 

where  v(k)  is  the  ra-diraenslonal  zero-mean  white  Gaussian  discrete- 
time  noise  with 

E{v(i  )vT( j  )  }  -  RS  t j  (24) 

and  assumed  independent  of  both  x(0)  and  w(*). 

The  system  of  Equations  21-24  describes  a  time-invariant  system 
with  stationary  noise.  If  the  system  model  is  stabillzable  and 
detectable,  there  will  be  constant  gain  for  both  the  controller  and 
observer  dynamics.  The  stochastic  optimal  control,  u,  will  minimize 
a  quadratic  cost  function  of  the  form: 


J  -  E{Va  l  CxT(i)$x(i)  ♦  uT(i)£u(i)]}  (25) 

1-1 

The  LQG  solution  for  the  time- invar iant  system  with  the  station¬ 
ary  random  noise  model  is^11*  12^ 

u (k )  -  -G*x (k )  (26) 

g*  -  (6+rTPr)-1rTP*  (27) 


where  the  P  satisfies  the  algebraic  Ricatti  equation 

P  -  ( 4-TG* )TP( *-TG*)  ♦  G*  TR  G*  +  3  (28) 

The  state  estimate,  x(k),  is  given  by 

x ( k )  -  x (k/k-1 )  +  K ( Y (k  )  -  H$(k/k-1))  (29) 

where 

x(k/k- 1 )  -  $S(k-l)  +  Tu(k-I)  (30) 

and 

K  -  PeHT(HPeHT  ♦  R)-1  (31  ) 

The  covariance  of  the  state  estimate,  PQ,  satisfies  the  algebraic 
Riccati  equation 

Pe  "  (*“K*H)Pe($-K*H)T  +  *KRKT«T  +  Q  (32) 


14 


where 


K*  -  OK  (33) 

2.4.2  The  Separation  Property  and  Suboptlmal  Design 

It  Is  wise  to  seek  the  possibility  of  a  separate  determination  of 
the  optimal  state  estimate  and  the  optimal  controller.  Due  to  com** 
putional  delay,  what  is  often  considered  for  implementation  is  not 
the  control  law  given  by  Equation  26,  but  a  suboptlmal  control 

u (k )  -  -G*x (k/k- 1 )  (34) 

wfrere 


x(k  +  1 /k) 


®x(k/k-1 )+ru(k)  +  0K(Y(k )-Hx(k/k-1 ) ) 


Define  the  error  of  the  predicted  estimate  of  the  state  as, 
Ax (k/k- 1 )  -  x (k )-x(k/k- 1 ) 


(35) 


the  x(k+1)  and  Ax(k+1/k)  may  be  combined  into  a  simple  system  de¬ 
scribed  by  the  state  equat i ons  I- 23  J 


'x  ( k ♦  1  ) 

^Ax(>*1  /k) 
where 


♦  -TG*  ToAlxik) 


■A 


«-K*H  /\x(k/k-  1  )/ 


(36) 


K’ 


♦  K 


The  output  equation  of  the  complete  system  has  the  form 


Y(k)  -  (H 


0) 


v  ( k  ) 


(37) 


From  Equation  36,  the  LQG  suboptimal  control  dynamics  also  consist 
of  two  parts:  the  dynamics  of  the  plant  with  a  feedback  controller 
and  the  dynamics  of  the  estimator  feedback  loop.  The  matrix  describ* 
ing  the  dynamics  of  the  closed  loop  controller  is 

A*  -  e-TG*  (38) 

whereas  the  analogous  matrix  for  the  closed-loop  estimator  is 

e-K»H  (39) 


A* 


where 


The  eigenvalues  of  A*  and  t*  are  independent  of  each  other.  In  the 
synthesis  of  the  control  loop  it  is,  therefore,  possible  to  arrange 
the  poles  of  the  estimators  and  controllers  separately.  Therefore, 
the  eigenvalues  of  the  closed-loop  system  simply  consist  of  the  regu 
lator  poles  and  observer  poles,  i.e.,  the  eigenvalues  of  A*  and  A*, 
respectively.  It  is  well  k-nown  that,  in  the  time- invari ant  system 
case,  if  and  only  if  the  open-loop  system  model  is  both  stabilizable 
with  respect  to  u(k)  and  detectable  wi^n  respect  to 

y(k),  there  do  exist  gains  G*  and  K*  tnat  can  provide  asymptotic 
closed-loop  stability.  Furthermore,  under  the  stronger  assumptions 
of  complete  controllability  and  complete  observability,  we  can  place 
both  the  regulator  and  observer  poles  arbitrarily  (within  the  re¬ 
striction  of  the  complex  conjugate  pairs). 

It  is  assumed  for  system  (21),  (23)  that  the  pair  ($,  r)  is  either 
completely  controllable  or  stabilizable  by  state  feedback,  and  the 
pair  ($,  H)  is  either  completely  observable  or  detectable,  then  the 
gain  matrix  of  the  controller,  G*  ,  and  the  gain  matrix  of  observer, 

K* ,  may  be  obtained  from  Equations  27,  28  and  31,  32,  33,  respective 
ly.  They  are  functions  of  the  weighting  matrices  R,  §  and  the  obser 
vational  noise  variance,  R,  and  the  system  noise  variance,  Q,  respec 
tively.  Therefore,  the  eigenvalues  of  the  closed-loop  controller  A* 
will  be  changed  by  means  of  the  changes  of  the  weighting  matrices, 

R,  and,  similarily,  the  eigenvalues  of  the  closed-loop  observer 
S*  are  also  changed  with  the  variation  of  the  matrices,  R  and  Q. 

Let 


Q  *=  HtH  ,  R  -  uoI0 
Q  ■  uqIo  ,  R  -  urIo 

In  order  to  find  the  appropriate  arrangement  of  the  controller  and 
observer  poles,  it  is  necessary  to  study  the  loci  of  the  A*  eigen¬ 
values  and  the  A*  eigenvalues  with  uq  or  y  r  ,  sjq  .  The  loci  of  the 
maximum  and  minimum  moduli  of  the  eigenvalues  of  A*  for  the  systems 
of  Case  6  and  Case  7  in  Table  2  are  listed  in  Table  3. 

A  x 

The  maximum  and  minimum  moduli  of  the  eigenvalues  for  A  vs.  the 
different  parameters,  ur  and  yq ,  are  listed  in  Table  4. 

It  is  assumed  that  the  covariances  of  the  measurement  noise  and  the 
system  (plant)  noise  are  Rg  -  1 0  ^  I  q  and  Qg  -  1 0  “* 1  0  *  respectively. 
Based  on  this  data  and  Tables  2-4,  we  can  arrange  the  position  of 
the  controller  and  observer  poles.  The  procedures  are  as  follows: 

(1)  Placement  of  actuators:  Based  on  the  calculation  of  the  degree 
of  cont r ol lab  i  1 i t y  for  different  arrangement  of  actuators,  to  deter¬ 
mine  the  best  placement  of  the  actuators.  As  shown  in  Table  2,  if 
the  number  of  actuators  is  12,  the  better  placements  are  indicated 
in  Cases  7  and  8.  Case  7  is  proposed  to  be  simulated  for  the  design 
of  the  controller  and  observer. 


TABLE  3.  THE  RELATIONSHIP  BETWEEN  THE  POSITION  OF  THE  REGULATOR 

POLES  AND  WEIGHTING  MATRIX  R 


A 

R  -  uclo 
UC 

Case 

Moduli  of  A* 
Min. 

6. 

Eigenvalues 

Max  . 

Case  7 

Moduli  of  A* 
Min. 

Eigenvalues 

Max . 

10-16 

0 . 3071  7x10-7 

0.99997 

0.30334x10-7 

0.99998 

10“  1  * 

0.42290x10-7 

0.99969 

0.33982x10-7 

0.99982 

1  O'*1  2 

0.1  1  522x 1 0“5 

0 .99696 

0.39831X10*6 

0.99826 

10-10 

0.11217x10^3 

0.97817 

0.36832x10“u 

0.98271 

10“8 

0.1  1  T 1  8x  1  0  "* 1 

0.97818 

0. 3731 6x10-2 

0.83858 

1  O'*6 

0.48849 

0 . 97865 

0.35403 

0.93776 

10-4 

0.90319 

0 . 98828 

0.86382 

0.98745 

10‘2 

0.96867 

0 . 99848 

0.96763 

0.99806 

1 

0.98768 

0.99984 

0.98809 

0.99981 

1  02 

0.99037 

0.99998 

0.99038 

0.99997 

1  o1* 

0.99042 

0.99999 

0.99042 

0.99999 

1  o6 

0.99042 

0.99999 

0.99042 

0.99999 

TABLE  4.  THE  MAX.  AND  MIN.  MODULI  OF  THE  EIGENVALUES  OF  A«  VS.  PARAMETERS  uR,  Uq  FOR  CASE  6  AND  CASE  7 
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(2)  Placement  of  observer  poles:  Based  on  the  values  of  the 

measurement  noise  and  plant  noise  covariances,  i.e.,  Rg  -  Qs 

-  IO^OIq,  the  minimum  and  maximum  moduli  of  the  A*  eigenvalues  may 
be  determined  preliminarily  from  Table  4.  In  general,  as  for  the 
design  of  the  Kalman  filter,  the  R  and  Q  should  satisfy  the 
conditions:  R  >  Rg,  Q  >  Qg .  After  these  conditions  have  been 

satisfied,  the  location~of  the  observer  poles  may  be  changed  by 
means  of  the  variation  of  parameters  up  and  uq  in  Table  4.  When  R  = 
Rg,  Q  ■>  Qg ,  the  minimum  and  maximum  moduli  of  the  eigenvalues  of  A* 
(observer)  are  0.74399  and  0.99450,  respectively. 

(3)  Placement  of  controller  poles:  As  we  know^23]  the  minimum 
modulus  of  the  observer  closed^loop  eigenvalues  should  be  less  than 
that  of  controller  closed^loop  eigenvalues,  so  that  the  estimator 
can  provide  accurate  timely  state  information  for  the  controller. 
Based  on  this  principle,  the  possible  appropriate  modulus  of  the 
minimum  controller  eigenvalues  for  Case  7  (Table  3)  is  0.86382, 
i.e.,  R  «  10-1*lo,  the  minimum  and  maximum  modulus  of  the  eigenvalues 
of  A*  (controller)  are  0.86382  and  0.98745  (Table  3). 

(4)  Certifying  by  simulation:  In  order  to  certify  the  preliminary 
design  of  the  controller  and  observer,  in  adition  to  showing  the  LQG 
designed  controlled  response  for  this  system,  we  also  need  to 
propose  another  kind  of  controller  and  observer  placement  for 
comparing  with  the  above  preliminary  design. 

2.5.  Digital  Simulation 

The  simulation  of  the  LQG  digital  control  of  the  orbiting  shallow 
spherical  shell  system  is  considered  here  including  some  controllers 
and  observers.  The  software  package,  0RACLS,^24]  ha3  been  used  for 
implementation. 

It  is  assumed  that  the  measurement  accuracy  of  the  displacement 
sensors  is  about  the  order  of  1  centimeter  for  the  shell  deflection, 
and  the  modelling  error  for  the  dynamical  system  is  less  than  the 
error  of  the  measurement  sensors.  The  covariance  matrices  used  for 
the  simulation  of  both  the  sensor  measurement  and  system  (plant) 
noise  are  determined  as  follows: 

Rg  =  10-6iq  is  the  covariance  of  the  measurement  noise  used  in  the 
simulation  and  Qg  =  10“10Io  is  the  covariance  of  the  system  (plant) 
noise  used. 


Now  let  w(i)  be  the  system  noise 
then 

E ( w (  1  )  }  -  0  and  E{v(i)vT(j)} 

Define  v(i)  as  the  observational 
model .  Similarly, 

E { v ( 1  )  }  -  0  and  e{v(1 )vT(J  )  } 


used  in  the  Kalman  filter  model, 

-  QSij 

noise  used  in  the  Kalman  filter 

-  R5iJ 
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Initial  state  simulation  is  then 

x(0)  -  x(0)  +  wq  x(0)  «  e{x(o)} 

Po  -  e{w0w0t} 

The  initial  conditions  are  assumed  as 

x  t  ( 0 )  -  <J>(0)  -  0.0  ,  x2<0)  -  ♦(O)  -  0.1  ,  x  3  ( 0 )  -  0(0)  -  0.1 

xjj(0)  -  e  i  ( 0 )  •  0.01  ,  X5(0)  -  xg(0)  •>  ...  -  x^g(0)  ■  0.0 

and  the  initial  values  of  the  estimated  state  are  assumed  to  be 
zer o  ,  i  .  e  . , 

x  t ( 0 /- 1  )  -  0.0  ,  x2(0/-1)  -  0.0  , 


X3<  0/-1  )  *  ...  ■  x^g(0/-1)  ■  0.0. 

In  order  to  study  the  influence  of  different  numbers  of  actuators  on 
the  LQG  control  process  and  the  influence  of  different  locations  for 
a  given  number  of  actuators  on  the  LQG  control  process,  Cases  2,  4, 

6,  7  (Table  2)  are  selected  for  simulation.  The  minimum  modulus  of 
the  eigenvalues  of  A*  for  the  four  cases  are  0.91717,  0.86701, 
0.90319,  0.86383,  respectively,  but  the  minimum  modulus  of  the  eigen¬ 
values  of  A*  are  the  same,  i.e.,  0.74399.  The  transient  response 
under  the  LQG  optimal  digital  control  for  Cases  2,  4,  6,  7  are  shown 
in  Figures  2-5.  It  is  evident  that  the  transient  response  of  Case  7 
is  the  best  among  Cases  2,  4,  6,  7,  because  there  are  more  actuators 
for  Case  7  ajid  the  degree  of  controllability  for  Case  7  is  also  the 
highest  among  the  four  cases.  It  is  valuable  to  point  out  the  fact 
that  the  transient  response  of  Case  6  is  worse  than  that  of  Case  4 
even  though  the  number  of  actuators  in  Case  6  is  greater  than  those 
in  Case  4.  The  reason  is  that  eight  actuators  for  Case  6  are  all 
located  on  the  same  nodal  circle  (£  »  const.)  and  the  influence  of 
the  control  forces  along  the  same  nodal  circle  on  the  displacement 
at  the  other  nodal  circles  are  very  weak.  In  view  of  the  degree  of 
controllability,  it  is  also  easily  understandable.  The  degree  of 
controllability  of  Case  6  is  less  than  that  of  Case  4. 

In  summary,  the  locations  of  the  actuators  should  be  as  far  as  possi¬ 
ble  from  the  nodal  meridians,  and  also  should  be  arranged  so  that  as 
few  actuators  as  possible  will  be  located  on  the  same  nodal  circle. 

As  we  pointed  out  in  Ref.  [23],  an  appropriate  design  of  LQG  optimal 
control  systems  should  require  that  the  minimum  modulus  of  the  eigen¬ 
value  of  the  LQG  closed-loop  observer  should  be  less  than  the  mini¬ 
mum  modulus  of  the  eigenvalue  of  the  LQG  closed-loop  controller,  so 
that  the  estimator  can  provide  accurate  timely  state  information  for 
the  controller.  In  order  to  certify  again  this  fact,  four  combina¬ 
tions  of  observer/controller  pole  locations  (EC-1,  EC-2,  EC-3,  EC-4) 
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Figure  2.  LQG  Control  of  Shallow  Spherical  Shell  Transient  Response  for  Case 
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LQG  Control  of  Shallow  Spherical  Shell  Transient  Response  for  Case 
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Figure  5.  LQG  Control  of  Shallow  Spherical  Shell  Transient  Response  for  Case 


described  in  Table  5  are  proposed  for  simulation  based  on  Case  7  (12 
actuators,  see  Table  2). 

The  transient  responses  of  the  LQG  optimal  digital  control  of  the 
shallow  spherical  shell  for  EC-1  to  EC-il  are  shown  in  Figures  6*10. 

It  is  evident  that  the  transient  response  of  attitude  and  modal 
amplitude  for  EC-1  are  the  worst  (Figures  6  and  7).  These  transient 
responses  do  not  converge,  because  the  minimum  modulus  of  the 
observer  eigenvalues  is  greater  than  the  minimum  modulus  of  the 
controller  eigenvalue;  thus  the  estimator  cannot  provide  the  timely 
accurate  estimate  of  the  state.  The  transient  response  of  EC- 2  is 
better,  because  the  minimum  modulus  of  the  observer  eigenvalues  Is 
less  than  that  of  the  controller,  so  the  observer  can  provide  the 
accurate  timely  state  estimate  (Figure  8).  It  is  expected  that  EC-3 
and  EC-4  both  can  provide  appropriate  timely  state  estimates  for  the 
controller,  (Figures  9,  10)  but  the  terminal  part  of  the  transient 
response  of  the  estimate  and  control  for  EC-r  is  degraded  (note  the 
lack  of  smoothness  in  Figure  10).  The  appearance  of  this  phenomenon 
is  due  to  the  small  ratio,  ur/uq,  so  that  Kalman  filter  becomes  too 
sensitive  to  the  new  observational  data  with  the  actual  measurement 
noise.  Therefore,  the  relationship  between  the  locations  of  the 
controller  and  observer  poles  for  EC-2  and  EC-3  are  better. 

2.6.  Conclusions 

The  analysis  and  design  of  the  optimal  LQG  digital  shape  and  orienta- 
tion  control  for  an  orbiting  shallow  spherical  shell  system  are  pre¬ 
sented.  The  emphasis  in  this  paper  is  placed  on  the  analysis  and 
design  of  LQG  optimal  digital  controllers  and  observers  for  the 
orbiting  flexible  shallow  spherical  shell  system.  As  for  the  place¬ 
ment  of  the  controller  and  observer  poles  the  minimum  modulus  of  the 
eigenvalue  of  the  closed-loop  observer  must  be  less  than  the  minimum 
modulus  of  the  eigenvalues  of  the  closed-loop  controller,  so  that 
the  observer  can  provide  the  accurate  estimate  of  the  state 
variables  for  the  controller.  When  placing  the  positions  of  the 
controller  and  observer  poles  to  meet  the  above  requirement,  we 
should  also  pay  attention  that  the  ratio  ur/uq  cannot  be  too  small; 
otherwise  the  Kalman  filter  becomes  too  sensitive  to  the  new 
observational  data  so  that  state  estimate  will  not  be  so  smooth. 

The  problems  of  determining  the  number  and  location  of  the  actuators 
are  also  studied  by  means  of  the  concept  of  the  degree  of  controlla¬ 
bility  together  with  transient  response  simulations.  For  the  model 
of  the  orbiting  shell  treated  here  12  actuators  properly  placed  can 
provide  both  shape  and  orientation  control.  Actuators  placed  along 
the  edge  of  the  shell  should  be  capable  of  providing  thrust  both 
tangentially  and  normal  to  the  base  of  the  shell. 
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TABLE  5.  THE  PARAMETERS  OF  THE  SYSTEM  EC1  -  EC* 


Combination  of 

yv 

Min.  Modulus 

Min.  Modulus 

Obser./Contr. 

R  -  UCI0 

R  -  prIq 

Q  -  WqIo 

of  Eigenvalues 

of  Eigenvalues 

Poles 

PC 

PR 

PQ 

of  A* 

of  ft* 

EC-1 

10-4 

10“6 

10-16 

0.86383 

0.99*02 

EC-2 

10-4 

10-6 

10**10 

0.86383 

0.7*399 

EC-3 

10-4 

10**6 

10“6 

0.86383 

0.32581 

EC-* 

IQ'4 

10“6 

10-4 

0.86383 

0.1*8*7x10“2 
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Figure  6.  LQG  Control  of  Shallow  Spherical  Shell  Transient  Response  for  System  EC-1 (a) 
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fG  Control  of  Shallow  Spherical  Shell  Transient  Response  for  System  EC-2 


Attitude  Angle  (rad.) 
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Figure  9.  LQG  Control  of  Shallow  Spherical  Shell  Transient  Response  for  System  EC-3 
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Figure  10.  LQG  Control  of  Shallow  Spherical  Shell  Transient  Response  for  System  EC-4 


3.  SOME  DEFINITIONS  OF  DEGREE  OF  CONTROLLABILITY  (OBSERVABILITY) 

FOR  DISCRETE-TIME  SYSTEMS  AND  THEIR  APPLICATIONS 

3.1.  Introduction 

The  concept  of  controllability  is  one  of  the  cornerstones  of  modern 
control  theory.  Yet  despite  its  fundamental  importance  from  a  theo¬ 
retical  point  of  view,  its  practical  utility  for  control  system 
either  controllable  or  it  is  uncontrollable.  There  is  no  provision 
for  consideration  of  the  more  subtle  question:  how  controllable  is 
the  system? 

A 

The  desirability  of  a  degree  of  controllability  (earlier  it  was 
called  the  controllability  index)  concept  has  been  recognized  in  the 
literature  since  1961  when  Kalman,  Ho  and  Narendra^2^]  discussed  it. 

Early  papers  in  the  area^2^“2®^  concentrated  on  particular  proper¬ 
ties  of  either  the  controllability  Grammian  matrix  or  the  controlla¬ 
bility  matrix  itself  in  developing  definitions  of  the  degree  of 
controllability.  It  is  natural  to  try  to  connect  the  degree  of  con¬ 
trollability  to  properties  of  the  standard  controllability  matrix 
Pc  =  ( B :  A  B  :  ...  An* ’ B) ,  and  define  the  degree  of  controllability  as 
the  square  root  of  the  minimum  eigenvalue  of  PCPCT.  In  1979  Viswana- 
than,  Longman  and  Likins  considered  the  "shortcomings"  oi  is  defi¬ 
nition^2^-!.  Four  apparent  difficulties  with  this  definition  must 
somehow  be  handled  before  the  definition  becomes  "viable";  there¬ 
fore,  they  presented  a  new  definition  of  the  degree  of  controlla¬ 
bility  of  a  control  system  by  means  of  a  scalar  measure  based  on  the 
concept  of  the  "recovery  region." 

In  this  paper  three  candidate  definitions  of  degree  of  controlla¬ 
bility,  which  were  presented  first  by  Muller  and  Weber^2®^  for 
linear  continuous  systems,  are  presented  for  linear  discrete-time 
systems  based  on  the  scalar  measure  of  the  Grammian  matrix.  The 
three  candidates  for  the  degree  of  observability  of  linear  discrete¬ 
time  systems  are  also  presented.  Because  some  difficulties  with 
this  definition  (as  pointed  out  by  V i swanathan t 29 ]  have  been  handled 
here,  and  the  degree  of  controllability  (observability)  based  on  a 
scalar  measure  of  the  Grammian  matrix  can  be  readily  calculated,  the 
three  candidates  for  the  degree  of  controllability  (observability) 
are  viable  for  practical  engineering  design. 

The  emphasis  of  this  paper  is  in  showing  the  physical  and  geometri¬ 
cal  meanings  and  the  general  properties  of  the  three  candidate  defi¬ 
nitions  of  degree  of  controllability  (observability)  as  physically 
meaningful  measures.  The  concepts  of  degree  of  controllability  are  » 

applied  to  the  actuator  placement  problem  for  the  orbiting  shallow 
spherical  shell  control  system.  The  degree  of  controllability  for 
seven  cases  In  which  the  placement  of  actuators  are  all  different  * 

are  compared.  The  LQG  transient  responses  for  several  typical 
systems  with  different  degrees  of  controllability  are  also  shown. 

3.2.  Concept  and  Physical  Meaning  of  Controllability 

In  order  to  show  the  physical  meaning  of  the  degree  of  controlla¬ 
bility,  first  of  all,  we  consider  the  fuel  optimal  problem  for  the 
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linear  discrete-tine  system.  The  state  equations  are  as  follows: 

x ( k i )  -  •(ki,ki-1)x(k1.1)  +  r(k1-1 JuCki-! )  (40) 

(i  -  1 ,  2,  ...  N ) 

The  fuel  optimal  control  problem  can  be  stated  as  follows:  To  find 
the  system  control  ufk^j)  (i  -  1,  2,  ...  N)  which  can  transfer  from 
the  initial  state  x(ko)  to  the  final  state  of  the  system  x(kN)  -  0 
during  the  time  interval  [k0,kfj],  and  such  that  the  performance 
i  ndex , 


N 

JN  -  l  <  u ( k i - i ) ,u(ki-i )  >  (41) 

i  - 1 

is  minimized . 

Equation  (40)  may  be  written 

N 

x(kN)  -  «(kN,k0)x(ko)  +  l  ♦(kfl.kj)  T(ki-1)  u(kj-i)  (42) 

i-1 


or 


x(kN)  -  <Kku.ko)  x(ko)  -  TpU 


(43) 


where 


Tp  -  (  r(kk-i  )  :  ♦(kfj,  k  p|  _  i  )  T  (  kfj-2 )  : 

fu  (k(|-i  ' 
u(kM-2) 


U  - 


[u(k0)  J 


♦  ( k(j ,  k0)  f(k0) ) 


Considering  the  requirement  that  x(kfj)-0,  Equation  (43)  may  be 
wr i tten : 


(44) 


(45) 


-x(k0)  -  PCU  (46) 

where 

Pc  *  ®(ko»'<N)  Tp  (47) 

and  the  fuel  optimal  problem  becomes  the  following  conditional  extre¬ 
mum  problem 

inf  ||  |U||  2|  PCU  ( • )  -  -x(ko)}  (48) 

where  |  |  •  |  |  represents  the  Euclidian  norm.  This  problem  can  be 
solved  by  using  the  Lagrange  multiplier  method.  If  X  is  the 
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Lagrange  multiplier,  the  problem  stated  in  Equation  48  is  equivalent 
to  the  general  problem  of  extremizing  the  following  coat  function 


Jp-  <  U ,  U  >  -  XT( PcU  ♦  x(k0) ) 


In  order  for  Jp  to  be  an  extremum  the  following  necessary  conditions 
must  be  satisfied: 


f±»  -  2U-P0T*  .  0 

-  PcU+x (kg)  -  0 


(49) 


i  .e .  , 


(PCPCT) A  -  -2x(k0) 
2U  -  PcTA 


(50) 


If  the  system  is  controllable,  then  the  inverse  of  the  matrix, 
(PCPCT),  exists,  and  the  optimal  control,  U*,  i3 

U*  -  •*PcT(PcPcT)"*1x(ko)  (51) 


After  substituting  Equation  51  into  the  performance  index,  Equation 
41,  the  minimum  control  energy  is  obtained  as  follows: 

J*  -  <  U*.  0*  >  -  <  -PcT(PcPcT)-1x(ko)  ,  -pcT(pcpcT)-1x(ko)  > 


-  <  (PcPcT)'*1x(k0)  ,  x  ( k0 )  > 


(52) 


J*  -  xT(k0)Wc'*1  (kN,k0)x(kg) 


(53) 


where,  Wo(kN,k0)  -  PcpcT 

From  Intuitive  considerations,  the  linear  system  "1”  is  more  control¬ 
lable  than  the  linear  system  "2"  if  the  energy,  ,  required  by 

system  "1"  for  transferring  the  initial  state,  x(kg),  to  the  final 
state  x(k^)  -  0,  is  less  than  the  energy  required  by  the  second 
system,  J*2 ,  for  transfer  between  the  same  initial  and  final 

states  within  the  same  time  period,  or  stated  mathematically. 


xT  (  k  o  )  Wc  ]  (  kfj ,  kg  )  x  ( kg  )£xT  (  kg  )  WC2  (  k  ,  kg  )  x  ( kg  )  Vx(kg)  eRn  (54) 

- 1  - 1 

Because  Wci(kj|,kg)  and  Wc2(k(j,kg)  are  two  symmetrical  positive 
definite  matrices,  the  relationship  (54)  should  be  satisfied  for  any 
initial  condition,  x ( kg ) eR  n ;  this  implies  that 

Wc](kN,kg)  <  WC2(KN.kg) 


(55) 


or 


Wci(kfj,ko)  >  WC2(k(j,ko>  (56) 

Considering  the  symmetry  of  the  matrices,  Wci  and  WC2»  the  following 
relationships  should  be  satisfied: 


Xi(Wc])  <  Xi(Wc2) 

(i  - 

Xi(Wd)  >  X1(wc2) 

(1  ■ 

n ) 

(57) 

n ) 

(58) 

where  Xi(Wc2),  Xi(Wc1),  X  (Wc2 )  are  the  1th  eigenvalues  of 

1  1 

the  matrices,  Wci ,WC2»  Wc1  and  WC2«  respectively.  It  is  evident 
from  Equation  53  that  the  loci  of  equicontrol  effort  (J*  -  con* 
stant  is  a  superell i psoi d  and  its  equation  (J*  -  17  in  the  principal 
axis  system  is  as  follows: 


JLl 


i _ 12} 


+  ...  + 


¥  2 

_ 


1 


1  1 

xTTwT^T  x2(wcJ-< )  xn(wc* 1 ) 

Considering  the  relationship  between  Xi(Wc‘‘1)  and  Xi(Wc) 


Al(Wc"1}  “  Xn+1li(Wc) 


(59) 


(60) 


i  .e . , 


X^H^1)  <  x2(wc“1)  <  ...  <  Xn(Wc-1) 

A1(WC}  -  Xndc'1')  <  A2(WC>  -  Xn-1|wc'M')' 
<  .  .  .  <  Xn(Wc)  -  ^  (icJ-  I  ) 


Then  Equation  59  can  be  written  as  follows: 
X~2 


X -  2  x  2 

X! (wc)  +  x2(wc)  • • •  xn?wc) 


1 


(61  ) 


Therefore,  from  the  geometrical  viewpoint,  it  is  clear  from  Equa* 
tions  54  and  61  that  if  system  "1"  is  more  controllable  than  system 
"2,"  each  axis  of  the  superel 1 i psoi d  for  system  "1"  is  longer  than 
the  corresponding  axis  for  the  system  "2."  The  square  root  of  each 
of  the  eigenvalues  of  the  matrix,  Wc,  is  just  the  length  of  the 
corresponding  super ell i psoi d  principal  axis. 

3.3  Three  Candidates  for  the  Definitions  of  Degree  of  Controlla* 
bility  and  Their  General  Properties 

As  we  see  in  Equation  54,  in  terms  of  the  matrices  this  comparison 
of  the  costs  of  two  systems  leads  to  (56),  i.e., 
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WclO<N,ko)  >  WC2(kjj.ko) 


(56) 


which  means  that  Wci-Wc2  has  to  be  a  positive  semidefinite  matrix. 

The  matrix  condition  (56)  is  equivalent  to  "n"  scalar  conditions 
and,  therefore,  a  simple  scalar  quantity  cannot  be  employed  to 
describe  the  conditions  (56).  For  engineering  purposes  it  is  desir- 
able  to  replace  the  matrix  measure  by  a  scalar  figure  of  merit  [25]. 
There  are  three  obvious  candidates  for  such  scalar  measures: 
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the  maximum  eigenvalue  of  Wc  (kfl.kg),  the  trace  of  W0  (kfl,kg),  and 
the  determinant  of  Wc(kfj,ko)  [28].  The  significance  of  these  quanti- 
ties  is  shown  in  the  following  section: 

(1)  The  first  candidate  for  the  degree  of  controllability,  pi .  As 
shown  in  (54),  we  must  consider  the  values  of  the  control  effort 

xT(k0)WG  (kN,kg)x(kg)  Vx(kg)eRn  when  comparing  the  control  effort 
of  the  two  systems.  It  is  natural  to  use  the  maximum  control  ef - 
f  or  t , 

Max  xT(k0)Wc“1 (kN.ko)x(ko) 

|  |  x ( k o )  ||=1  (62) 

x(k0)eRn 

as  a  scalar  measure  of  the  control  effort.  It  is  well  known  that 
the  scalar  quantity  defined  in  (62)  is  just  the  maximum  eigenvalue, 

i  *  1  « i 

Amax(wc  )i  of  the  matrix  WG  ,  and  the  Amax(wc  )  13  Just  the  recipro¬ 
cal  of  the  minimum  eigenvalue  of  the  matrix,  Wc,  i.e., 

Max  xT(k0)Wc1 (kN,kg)x(kg)  ■  Amax(Wc  ’’  (63> 

I  I  x  ( k  o )  |  |  -  1 

x(k0)eRi 

where  the  Amin(wc^  is  the  minimum  eigenvalue  of  the  matrix, 

Wc  ( kjy ,  kg ) .  For  the  reasons  given  above,  we  may  define  the  degree  of 
controllability,  pi ,  as  follows: 

U1  ”  Ami  n  ^c  t  »  ^0  ))  “  °min^^c^  (64) 

This  definition  means  that  as  p^  increases,  then  the  control  effort 
decreases.  It  is  evident  that  the  square  root  of  the  degree  of  con¬ 
trollability  defined  by  Equation  (64)  is  just  the  minimum  semiaxis 
of  the  superell i psoi d  defined  by  Equation  (61)  (it  can  be  called  the 
degree  of  controllability  super  el 1 i psoi d  ) .  If  p^  =  0,  the  system 
will  not  be  controllable. 

(2)  The  second  candidate  for  the  degree  of  controllability,  p^. 
Similar  to  (62),  we  can  also  use  the  average  value  of  the  control 
effort,  xT(kg)Wc'‘1x(ko)  ,  over  the  unit  hypersphere  fx  (kg)  :|  |  x  (kg)|  | 

-  1  } 
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xT(ko)Wc'‘1x(ko)dx(ko) 

x  ( k  n  )1  I  -  1  _ 


I  I  x  ( k  o  )|  | 


dx(k0) 


1 


as  a  scalar  measure  of  the  control  effort.  After  Integrating,  it 
can  be  obtained  [28] 


■  M  x  (k0)|  |  -  1 

dx(k0) 

|  |  x  (kO  )|  |  -  1 


xT(k0)Wc“1x(k0)dx(ko) 


-trV  *1 
ntrwc 


(65) 


where  n  is  the  order  of  the  state.  Therefore,  we  may  define  the 
second  candidate  for  the  degree  of  controllability  as  follows: 


U2  “  trWc-'(kN,k0) 


(66) 


For  practical  applications  it  is  desirable  to  maintain  the  average 
cost  Equation  (65)  as  small  as  possible.  Hence,  the  measure,  V2> 
has  to  be  as  large  as  possible.  The  definition  (66)  of  a  measure 
instead  of  (65)  is  more  convenient  because  uncontrollable  systems 
are  characterized  by  a  vanishing  value  of  U2>  which  arises  from  a 
limiting  process. 

(3)  The  third  candidate  of  degree  of  controllability,  m.  The 
third  possible  candidate  for  a  scalar  quantitative  measure  of  con¬ 
trollability  is  the  determinant  of  Wc(k«.ko)  because  the  volume  of 
controllability  superel  1  i  psoid  xT(ko)Wc  ^  (kyj  ,ko  )x  (kg)  *  1  is  propor- 
tional  to  the  square  root  of  the  determinant  of  Wc(kN,ko): 

Vol  =  [  dx (kO ) =Const .  /detWc  (67) 

x(ko)TWc"‘1x(ko)  <  1 

Therefore,  we  may  define  the  third  candidate  definition  of  degree  of 
controllability  as  follows: 

u3  =  (detWc(kN,k0) ) 1 /n  (68) 

From  geometrical  considerations,  2n(u3)n^2  is  just  the  volume  of  the 
hy per r ectangul ar  parallepiped  whose  sides  are  the  axes  of  the  hyper-1 
ellipsoid.  The  larger,  y3,  the  more  controllable  the  system  is;  the 
uncontrollable  system  is  also  characterized  by  u3  -  0. 

(*<)  The  general  properties  of  the  degree  of  controllability.  In 
the  previous  sections  we  gave  the  definitions  of  degree  of  controlla¬ 
bility  and  their  physical  and  geometry  meaning.  The  degrees  of  con¬ 
trollability,  as  a  scalar  measure,  are  invariant  under  an  orthogonal 
similarity  transformation.  The  proof  will  be  given  in  the  Appendix 
of  this  paper.  In  addition,  they  are  scalar  quantities  defined  over 
the  cone  set,  Wc  ,  of  positive  definite  and  semidefinite  matrices. 
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Their  properties  are  as  follows  [28]: 

(a)  Mi(Wc)  -  0  (1  -  1,  2,  3)  for  WceWc*  with  det  Wc  -  0,  l.e.,  the 
system  cannot  be  controllable; 

(b)  m (Wc )  >  0  (1  -  1,  2,  ...  N )  for  WceWc*  with  det  Wc  >  0; 

(c)  y i ( A Wc )  -  Xyi(Wc)  (1  -  1,  2,  3)  for  X  >  0,  WceWc*; 

(d)  mi(Wc)  >  u i  ( Vfc -|  )  +  Ui(WC2)  (i  ■  1,  2,  3)  for  We  ■  Wci  ♦  WC2» 
W01.  WC2eWc*. 

The  proof  of  the  properties  (a)  -  (d)  will  be  given  in  the  Appendix 
of  this  paper. 

The  physical  meaning  of  the  properties  (a)  -  (c)  are  evident.  The 
physical  meaning  of  property  (d)  is  that  the  degree  of  controlla- 
bility  for  the  system  which  contains  two  control  elements  Is  higher 
than  the  sum  of  the  degrees  of  controllability  for  the  two  systems 
each  of  which  contain  only  a  single  controller. 

The  ui .  U2*  P3  are  three  kinds  of  candidates  for  the  degree  of  con¬ 
trollability,  the  physical  meanings  of  which  are  different.  The 
following  relationship  among  u-|  ,  U2  and  ^3  can  be  stated  [28]: 

ui  <  U2  i  U3  (69) 

The  sufficient  and  necessary  condition  for  the  equality  to  be  true 
Is  that  all  of  the  eigenvalues  of  Wc  be  equal. 

3.1*.  The  Definition  of  Degree  of  Observability  and  Its  Physical 
Meaning 

Suppose  the  dynamical  equations  and  observational  equations  are  as 
f ol 1 ows : 

x(kt)  -  $(ki ,ki»i )x  ( k i ^  1 )  +  r(k1*i )u(ki-i )  (70) 

y  ( k  i-.  1  )  -  H(ki«1 )  x  ( k  i-  1 )  ♦  v(ki-i)  (1  -  1,  2,  ...  N)  (71) 

where  v(k^l)  Is  observational  white  noise  with  statistical  proper¬ 
ties 

Ev(kj«,})  -  0  E  { v  ( k  1  1  )  vT  (  k  j  j  (1  -  1,2,  ...  N). 

Equation  (70)  may  be  written  in  another  form  in  terms  of  the  initial 
state  as  follows: 


i 

x^i)  -  ♦<kliko)x(k0)  ♦  l  ♦(kilkj)r(kj-1)u(kj-1) 

y(kA<1)  -  H(k^i  )*(ki*1  ,k0)x(k0) 


(72) 
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(73) 


i-1 

♦  HCki-l)  l  ♦(k1-i ,kj ) T ( k j - 1 )u(kj-! )  ♦  v(kj-^) 

J-1 

(t  -  1,2,  ...  N ) 

When  l  assumes  values  from  1  through  N,  Equation  (73)  can  be  put  in 
a  matrix  form  as  follows: 


Y  -  Pox(k0)  ♦  V 


(74) 


where 


Po  - 


f  y  (vc0  > 

y ( k i )  -  H(k1 ) r(ko)u(kg) 

N-1 

^y(kjj-«i)”*H(k}j-i)  £  $(k[j-i,kj)  r(kj.^)u(kj.^) 

/,H(k0)  \ 

H( k  i )  ♦(  k  ■,  ,k0) 

H(k{j-i  )*(kjj-i  ,ko) 


fv(k0)  ) 
v  ( k  i  ) 

V(kN-t)J 

E { V }  -  0  E{VVT} 


Ro 


Ri 


V 


rn-i 


We  consider  the  weighted  least  square  estimate  problem,  the  perform¬ 
ance  index  of  the  weighted  least  square  estimate  is 

J  -  (Y-Pox(k0))TR-1(Y-Pox(k0))  (75) 

It  is  well  known  that  the  weighted  least  square  estimate  x(ko)  is 

x(k0)  -  (P0TR‘1P0)-1poTR“17  (76) 

The  estimate  error  is 

x(k0)  -  x(k0)  -  ( P0tR“1 P0)-1 poR-1  (P0x(k0)-Y) 

-  -(P0TR“'IP0)“1P0,rR,*1V 
The  covariance  of  the  estimate  error  is 

E  { (  x  ( k0 )  -  x(k0) )  (x(k0)-x(k0)  )T}  -  (P0TR'*1Po^1  (77) 

I,  the  covariance  of  the  least  square  estimate  error  is 


When  R 
just 
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E  {  (  X  (kg )  -  x(kg)) (x(kg)*x(kg))T}  -  (P0Tpo^'*1  *  W g “ 1  (78) 

where 

w0  -  poTpo  (?9> 

Similar  to  the  degree  of  controllability,  we  call  the  Wg  the  degree- 
of  the  observability  matrix.  The  three  kinds  of  definitions  of 
degree  of  observability  are  as  follows: 

ui (Wg)  =  Amin(Wg)  (80) 

2  ( ^0 )  “  tpW*  I  <81  > 

U3(Wg)  «=  (det  Wg)  1  /n  (82) 


Because  the  Wg-1  is  the  covariance  matrix  of  the  least  square  esti* 
mate,  then  the  trWg1*1  is  the  sum  of  the  squares  of  the  variances 
associated  with  the  least  square  error  estimate,  so  the  physical 
meaning  of  y2(Wg)  is  evident.  ui(Wg)  is  the  maximum  eigenvalue  of 
the  covariance  error  matrix  Wg-*  for  the  least  square  estimte. 
u 3 ( Wg )  is  the  reciprocal  of  the  geometrical  average  of  the  error 
covariance  matrix  eigenvalues  for  the  least  square  estimate. 

In  general,  if  the  estimation  accuracy  of  system  "1"  is  higher  than 
that  of  system  "2,"  it  means 

WoT1  <  W02-1  (83) 

It  implies 

W01  >  W02.  trW01-1  <  trW02< 1  (84) 

and 

Xi(W01)  >  1 i ( W02 )  (i  =  1,  2,  ...  n)  (85) 

Therefore 

PI (Wqi )  >  ui ( Wg2 ) 

W2(Wg1  )  >  V*2  (  Wg2  ) 

U3 ( Wgi )  >  U3 ( Wg2 ) 

i.e.,  the  degrees  of  observability  defined  by  (80*82)  for  system  "1" 
are  all  higher  than  the  corresponding  degrees  of  observability  for 
system  "2." 

3.5.  The  Degree  of  Controllability  (Observability)  for  Discrete* 
Time  Invariant  Systems 

Let  the  dynamical  system  and  observational  system  for  the  discrete 
t i me* i nvar i ant  system  be  described  as  follows: 

x  ( k  i )  *=  *x(ki«i)  +  ru(kj[^l)  (86) 
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(87) 


% 


.  y  (k  i  )  -  Hx(kt)  +  v(kj[) 

The  degree  of  controllability  matrix  Wc  for  the  system  (86),  (87)  is 
Wc  -  PCPCT  (88) 


where 


Pc  -  (  T  :  ♦  T  :  ...  ♦n‘1D  (89) 

The  degree  of  observability  matrix  for  system  (86),  (87)  is 

Wo  “  P0Tpo  (90 ) 


where 


(91  ) 


The  three  kinds  of  definitions  of  degree  of  controllability  are  as 
follows : 


HI <wc )  -  »min<W0)  (92) 

“S'o’  •  tTJTT  193) 

U3(WC)  -  (det  tac ) 1 /n  (94) 

Similarly,  the  three  kinds  of  degree  of  observability  are  as 
follows : 

HI ( W0 )  -  *min(W0)  (95) 

y2(w0)  -  (96) 

U3(w0)  -  (det  W0)1/n  (97) 

3.6.  Application 


We  would  like  to  use  the  concepts  of  the  degree  of  controllability 
for  the  actuator  placement  problem  of  the  shallow  spherical  shell. 
The  motion  equations  for  a  shallow  spherical  shell  in  orbit  subject 
to  external  forces  and  torques  are  as  follows  [19]: 


X  -  AX  ♦  Bu  (98) 

where 


X  -  (7,  ♦,  0,  e i ,  e 2  •••  *  eg,  T,  *,  0,  ei,  e 2  ...  ,  eg)T 
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'iz2fy2+*y2fz2 
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lt  f  +  4  f  ^ 
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Jx*  ^U>Lic^ 
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yaw,  roll  and  pitch  angles,  respectively,  of  the 
undeformed  axis  of  the  shallow  spherical  shell. 

modal  amplitude  of  the  i^11  generic  mode  whose  shape 
function  is  $x  (  *  ) 

orbital  angular  rate,  constant  for  assumed  circular 
orbit 


4  charac ter i st i c  length  (the  base  radius) 

M i  i th  modal  mass 

<|>x(i)  the  x=axls  component  of  the  ith  modal  shape  function 

£xi»  *yi»  ^zi  the  components  of  the  moment  arm  for  the  ith  actuator 
Equation  98  is  nondimensional ized  according  to 
t  =  wct,  ej  =  qj(t)/A  (i  *=  1,  2,  ...,  6) 

The  derivative  in  Equation  98  is  with  respect  to  t. 


If  we  use  six  actuators  to  control  the  orientation  and  shape  of  the 
shallow  spherical  shell,  the  arrangement  of  the  six  actuators  may  be 
assumed  for  the  following  seven  cases  (Table  6): 
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TABLE  6.  THE  ACTUATOR  LOCATIONS  U,  8)  AND  FORCE  DIRECTIONS 

(fX,  f  y  »  fZ) 

CASE  1 

Actuator  No.  S  8  fx  fy  fz  Location  of  Actuator 

1  0.28 

2  0.57 

3  0.84 

4  0.28 

5  0.57 

6  0.84 


CASE  2 


Actuator  No. 

5 

8 

fX 

fy 

f 

1 

0.28 

5° 

1 

0 

0 

2 

0.57 

5° 

1 

i 

0 

3 

0.84 

5° 

1 

i 

1 

4 

0.28 

1  85° 

1 

0 

0 

5 

0.57 

1  85° 

1 

i 

0 

6 

0.84 

1  85° 

1 

i 

1 

Location  of  Actuator 


CASE  3 

Actuator  No. 

5 

8 

fx 

fy 

f  2 

Location  of  Actuator 

1 

0.28 

45° 

1 

0 

0 

2 

0.57 

45° 

1 

i 

0 

3 

0.84 

45° 

1 

i 

1 

4 

0.28 

225° 

1 

0 

0 

f — —  — j 

5 

0.57 

225° 

1 

1 

0 

V  y 

6 

0.84 

225° 

1 

1 

1 

CASE  4 


Actuator  No. 

5 

8 

f 

1 

0.28 

90° 

1 

2 

0.57 

90° 

1 

3 

0.84 

90° 

1 

4 

0.28 

270° 

1 

5 

0.57 

270° 

1 

6 

0.84 

270° 

1 

fy  fz  Location  of  Actuator 

0  0 
1  0 
1  1 
0  0 
1  0 
1  1 


CASE  5 

Actuator  No. 

5 

8 

fx 

f  y 

fz  Location 

i  of  Actuator 

1 

0.28 

45° 

1 

0 

0  _ 

2 

0.57 

45° 

1 

0 

0  xr^ 

r\ 

3 

1  .00 

45° 

1 

-sin45° 

cos45°  (  Vv 

_  ^ 

4 

0.28 

225° 

1 

0 

0  V 

V  ) 

5 

0.57 

225° 

1 

0 

0 

6 

1  .  00 

225° 

1 

s 1 n45  0 

-cos458  — 

— — 
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Table  6  (continued) 


CASE  6 


Actuator 

No. 

5 

8 

fx 

fy 

f  z 

1 

0.57 

0° 

1 

0 

0 

2 

0.57 

90° 

1 

0 

0 

3 

0.57 

1  80° 

1 

0 

0 

4 

0.57 

270° 

1 

0 

0 

5 

1  .00 

45  0 

1 

-sin45° 

cos  45 

6 

1  .00 

225° 

1 

Sin45 

0  - 

cos  45 

CASE  7 

Actuator 

No. 

5 

8 

f  y 

f  z 

1 

0.28 

90° 

1 

0 

0 

2 

0.57 

90° 

1 

0 

0 

3 

1  .00 

90° 

1 

- 1 

0 

4 

0.28 

270° 

1 

0 

0 

5 

0.57 

270° 

1 

0 

0 

6 

1  .00 

270° 

1 

1 

0 

Location  of  Actuator 


Location  of  Actuator 
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The  seven  cases  may  be  divided  into  two  groups.  The  first  group 
Includes  case  1  through  case  4.  The  second  group  includes  case  5 
through  case  7.  In  the  first  group  there  are  3ix  jets  (actuators) 
for  each  case.  In  the  second  group  there  are  six  jets  for  each 
case,  but  two  of  the  jets  for  the  second  group  are  located  at  the 
edge  of  the  shell  with  their  thrust  direction  tangent  to  the  shell's 
(circular)  edge.  In  order  to  reduce  the  possibility  that  the  jets 
used  primarily  for  shape  control  would  also  disturb  the  orientation 
of  the  shell,  the  placement  of  these  jets  are  arranged  symmetrically 
with  respect  to  the  shell's  undeformed  principal  axes. 

The  values  of  the  degree  of  controllability  for  each  case  are  listed 
in  Table  7.  The  reason  why  the  degree  of  controllability  for  case  1 
is  zero  is  that  all  the  actuators  are  located  along  the  meridional 
nodal  line  of  one  of  the  fundamental  shell  vibrational  modes. 

In  the  first  group  of  actuators,  the  best  placement  of  actuators  for 
which  the  degree  of  controllabililty  is  highest  is  case  4.  In  the 
second  group,  the  best  placement  of  actuators  is  case  7,  and  the 
degree  of  controllability  for  case  7  is  also  higher  than  that  for 
case  4.  The  reason  why  the  degree  of  controllability  for  case  6  is 
so  low  is  that  four  actuators  are  all  located  on  the  same  nodal 
circle  (£  *  constant). 

In  summary,  the  locations  of  the  actuators  should  be  as  far  as  possi¬ 
ble  from  the  nodal  lines,  and  also  should  be  arranged  so  that  as  few 
actuators  as  possible  will  be  located  on  the  same  nodal  circle. 

In  order  to  control  the  orientation  of  the  shell  effectively,  a  com¬ 
bination  of  tangential  jets  along  the  edge  of  the  shell  together 
with  selected  Jets  normal  to  the  shell's  major  surface  is  recom¬ 
mended. 

If  the  number  of  actuators  Is  limited  to  six,  the  arrangement  of 
actuators  as  in  case  7  is  suggested. 

The  transient  responses  with  the  LQG  control  for  case  2  and  case  7 
are  selected  to  show  the  differences  in  the  attitude  and  the  first 
three  modal  amplitude  responses. 

The  initial  conditions  are  the  same  for  the  two  cases,  i.e.,  roll 
*(0)  «  0.1  rad.,  yaw  4'(0)  *=  0.0,  pitch  0(0)  ■  0.1  rad.,  the  initial 
conditions  for  the  first  six  modal  amplitudes  are  qi(0)  -  1  meter, 
q2(t)  *  q3(0)  ...  -  qg(0)  *  0.0,  respectively. 

The  transient  response  of  the  attitude  motion  for  case  2  is  shown  in 
Figure  11,  the  transient  responses  of  the  1st  -  3rd  modal  amplitudes 
for  case  2  are  shown  in  Figure  12.  The  transient  responses  for  the 
attitude  motion  and  the  modal  amplitudes  for  case  7  are  shown  in 
Figure  13  and  Figure  14,  respectively.  The  comparison  of  the  first 
modal  amplitude  responses  for  case  2  and  case  7  is  shown  in  Figure 
1  5  . 
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TABLE  7.  THE  DEGREE  OF  CONTROLLABILITY  FOR  THE  DIFFERENT  CASES 


0.13875x10"12 

0.17646X10-11 

0.67227x10"? 

0.17673X10"1 1 

0. 23696x1 0“10 

0.25045x10"^ 

0.17974X10"11 

0.26434x10"10 

0.31558x10"6 

0.14332x10*10 

0.12421x10*9 

0.70448x10*6 

0.48360x1  O'*1  2 

0.46915X10"11 

0.61384x10*6 

0. 1 4332x1 0"10 

0.14788x10*9 

0.10084x10*5 

It  is  shown  in  Figures  11*15  that  the  roll  is  coupled  strongly  with 
yaw,  but  the  coupling  of  roll,  yaw  with  pitch  is  very  weak.  The 
quality  of  the  transient  responses  of  case  2  is  much  worse  than  that 
of  case  7.  This  is  a  result  of  the  differences  in  the  degree  of 
controllability  for  case  2  and  case  7. 

3.7.  Conclusions 

Three  candidate  definitions  of  the  degree  of  controllability  and 
observability  are  presented  for  linear  discrete* time  systems  based 
on  the  scalar  measure  of  the  Grammian  matrix.  Their  general 
properties,  together  with  the  physical  and  geometrical 
interpretations  for  the  fuel  optimal  control  problem  are  shown  in 
detail.  The  advantages  of  these  kinds  of  definitions  for  the  degree 
of  controllability  (observability)  are  the  clarity  of  the  physical 
and  geometrical  interpretations  and  the  simplicity  of  the  resulting 
calculation.  Thus,  they  are  very  useful  for  practical  engineering 
design. 

The  transient  responses  for  several  typical  systems  with 
different  degrees  of  controllability  show  that  the  quality  of  the 
system  transient  response  depends  on  its  degree  of  controllabl il ity : 

transient  responses  of  the  system  with  a  higher  degree  of 
controllability  are  better  than  those  of  the  system  with  a  lower 
degree  of  controllability.  The  applications  of  the  concept  of  the 
degree  of  controllability  for  actuator  placement  of  the  orbiting 
shallow  spherical  shell  system  are  successfully  implemented  in  this 
paper  . 


lMt 
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Figure  12.  LQG  Control  of  Shallow  Spherical  Shell  Transient  Response  for  Case  2  -  b(six  actuators) 


roll 
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Figure  13.  LQG  Control  of  Shallow  Spherical  Shell  Transient  Response  for  Case  7  -  a(six  actuators) 
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Control  of  Shallow  Spherical  Shell  Transient  Response  for  Case  7  b(six  actuators) 


amplitude  (meter) 
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Amplitude  for  Case  2  and  Case 


4.  OPTIMAL  DIGITAL  CONTROL  FOR  FREE*FREE  ORBITING  PLATFORMS 


4.1.  Specific  Alma 

To  develop  practical  design  methods  using  LQR  digital  controllers 
and  estimators  for  a  third  flexible  orbiting  plate,  focusing  on  the 
analysis  of  the  closed  loop  dynamics  of  the  platform  subject  to  dis* 
Crete  time  input  data. 

4.2.  Justification 

«  Large  flexible  spacecraft  systems  have  been  proposed  for  future 

applications  in  widespread  communication,  electronic  orbital  based 
mail  systems,  and  as  collector  of  solar  energy  for  transmittal  to 

•  Earth*based  receiving  stations.  For  such  missions,  with  the  inher* 

ent  size  and  necessary  low  weight  to  area  ratio,  the  flexible  parts 
of  such  systems  must  be  treated  as  nonrigid.  To  satisfy  the 
requirements  of  these  proposed  missions,  both  orientation  and  shape 
of  the  orbiting  system  should  be  controllable. 

Often  the  optimal  control  laws  for  this  system  are  developed  under 
the  assumption  that  the  state  vector  is  observed  directly  or  the 
state  information  can  be  stimated  on  a  continuous  basis.  However, 
for  future  applications  the  observational  data  will  often  be  col* 
lected  on  a  sampled  basis,  creating  a  discrete*time  data  system. 

Then  the  amount  of  information  collected  may  be  reduced  and  the 
format  of  data  input  may  be  acquired  more  conveniently.  For  the 
case  to  be  treated  here  only  a  deterministic  system  will  be  ad* 
dressed,  i.e.,  no  random  noise  will  be  considered. 

Therefore,  it  appears  useful  and  timely  to  3tudy  the  control  problem 
of  large  orbiting  space  structural  systems  with  discrete*time  obser* 
vational  data.  It  Is  well  known  that  the  development  of  modern 
control  theory  and  technology  provides  a  strong  tool  for  solving 
this  kind  of  engineering  problem.  C30,  31]  The  LQR  regulator  tech* 
nique^30,  31]  that  strong  tool  for  synthesizing  linear  system 
control  laws.  It  can  provide  acceptable  control  performance  once 
the  state  and  penalty  matrices  are  properly  selected.  It  does  not 
restrict  the  number  of  actuators  to  be  equal  to  the  number  of 
degrees  of  freedom  in  the  system.  Although,  the  LQR  method  has  been 
developed  and  widely  applied,  it  is  still  not  an  easy  task  to  apply 
it  to  the  engineering  design  for  the  control  of  large  space  structur* 
al  systems,  especially  for  systems  with  sampled  data  input.  There 
are  still  many  specific  problems  to  be  investigated.  These  are  the 

j  aims  of  this  proposed  research. 

4.3  Methodology 

4 

By  using  a  modified  version  of  the  general  formulation  of  the 
dynamics  of  a  general  flexible  orbiting  body  formulated  by 
Santini,L32]  the  equations  of  motion  for  a  free*free  beam  in  orbit 
were  developed.  In  order  to  gain  Insight  Into  the  dynamics  of  such 
a  large  flexible  system,  the  equations  of  motion  for  a  free*free 
beam  In  orbit  were  studied. t1]  The  motion  of  the  generic  point  in 
the  body  was  described  as  the  combination  of  the  rigid  body  motion 
plus  a  superposition  of  the  elastic  modes. 
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Assuming  the  center  of  mass  follows  a  circular  orbit  and  the  pitch 
and  the  flexural  deformations  occur  only  within  the  orbital  plane, 
it  is  seen  that  the  pitch  motion  does  not  influence  the  elastic 
motion.  Also,  the  pitch  and  the  elastic  modes  are  decoupled  for 
large  values  of  the  square  of  the  ratio  of  the  structural  modal  fre* 
quency  to  the  orbital  angular  rate.  For  small  values  of  this  ratio 
the  elastic  motion  is  governed  by  Hill's  B^term  equation  which  can 
be  approximated  by  a  Mathleu  equation.  Using  a  Mathieu  stability 
chart,  the  resulting  stability  was  considered.  For  small  amplitude 
flexural  motion,  the  rigid  body  and  elastic  modes  are  modelled  to 
the  first  order,  thus  linearizing  the  equations  of  motion. 

Subsequently,  an  extension  was  made  to  the  f orementioned  formulation  • 

and  was  applied  to  a  thin,  flexible  orbiting  plate.  [33] 

The  ability  to  accurately  determine  the  frequencies  and  mode  shapes 
is  essential  for  the  analysis  and  control  of  large  orbiting 
structures.  For  an  aluminum  square  plate  four  different  frequency 
and  mode  shape  approximation  methods  were  analyzed:^] 

1)  The  approximate  frequencies  and  mode  shapes  of  a  rectangular 
plate  were  derived  from  a  formulation  by  War  bur  ton .  [  31* ,  35] 

2)  The  analytical  results  for  a  square  plate  were  calculated  from  a 
method  developed  by  Lemke.^36] 

3)  The  frequencies  and  mode  shapes  were  computed  using  a  finite 
element  program,  STRUDL ,  written  at  M.I.T.L37] 

M)  The  frequencies  and  mode  shapes  were  computed  using  GT*STRUDL, 
an  updated  version  of  STRUDL  written  at  Georgia  Tech.^38] 

It  was  found  that  GT^STRUDL  obtained  better  results  than  STRUDL  and 
produced  accurate  results  for  specific  finite  element  input  grid 
point  locations  (node),  whereas  the  Warburton  and  Lemke  methods 
could  only  afford  approximate  answers. 

The  attitude  and  shape  control  can  now  be  achieved  by  placing  point 

thrust  actuators  perpendicular  to  the  main  surface  and  the  edge  of 

the  plate.  The  placement  of  the  actuators  on  the  main  surface  help 

control  the  shape  deformation  and  the  torque  about  two  axes.  The 

placement  of  the  actuators  along  the  edge  of  the  plate  help  control 

the  torque  about  the  third  axis.  Their  effects  on  the  rigid  body 

and  elastic  modes  are  modeled  to  the  first  order. C 39 D  for  this 

investigation,  it  will  be  assumed  that  the  sensors  and  actuators  are 

completely  colocated  with  the  actuators  and  that  the  system  is  , 

completely  observable. 

Now  the  control  laws  for  this  system  may  be  applied  to  obtain  the  < 

optimal  control  feedback  gains  based  on  an  application  of  the  linear 
regulator  problem  for  a  discrete^time  data  system. ^30,  31]  The 
implementation  of  the  LQR  procedure  will  be  accomplished  by  using 
the  ORACLS  routines  .  C24 ] 
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The  theorem  about  choosing  the  length  of  the  sampling  period  in 
order  to  guarantee  controllability  will  be  applied  to  select  the 
sampling  period. ^23]  When  the  controllability  of  the  continuous 
system  is  guaranteed,  in  order  to  maintain  controllability  for  the 
discretized  system  the  eigenvalues  of  the  system  must  satisfy 
certain  conditions.  As  for  the  signal  reconstruction,  the  sampling 
period,  AT  should  be  as  small  as  possible,  but  if  the  sampling  time 
is  too  small,  the  computational  requirements  may  exceed  the  computer 
speed . 

Under  normal  operation,  the  onboard  computer  estimation  and  control 
«  must  be  finished  processing  all  the  input  data  during  one  sampling 

period,  AT,  i.e.,  the  prediction  of  the  state  variable  which  will  be 
used  for  the  controller  must  be  estimated  in  time  before  the  be-1 
•  ginning  of  the  next  sampling  sequence.  Therefore,  the  sampling 

period  should  be  more  than  the  minimum  computational  time  required 
by  the  onboard  microcomputer  for  the  simulation  of  each  step  in  the 
estimation  and  control  process.  The  choice  of  sampling  time  is  also 
constrained  by  the  performance  of  its  transient  response,  i.e.,  over* 
shoot  characteristics,  settling  time,  steady  state  RMS  errors,  etc. 

4.4  Expected  Results 

A  comparative  parametric  study  will  be  performed  for  a  different 
sampling  period,  AT;  the  final  determination  of  the  sampling  period 
will  depend  on  the  compromise  between  several  factors  including  the 
controllability,  restoration  of  the  d i ser ete* time  input  signal,  the 
limitation  of  hardware  functions,  and  the  quality  of  the  transient 
responses . 

Software  such  as  that  currently  available  in  the  ORACLS  LQR  package 
can  be  used  to  synthesize  the  control  laws  and  simulate  the  trans* 
lent  responses  and  required  control  efforts. 

4.5  Why 

As  we  all  know  there  have  been  several  proposals  for  large  space 
structures,  such  as  the  Space  Station,  large  antennas,  solar  panel 
arrays,  propulsion  devices,  etc.  Often  these  modern  spacecraft 
systems  can  be  modeled  as  a  square  or  rectangular  plate  —  the  base 
—  with  the  attached  appendages.  The  acquired  knowledge  from  this 
proposed  study  will  have  a  wide  range  of  applicability  in  the  area 
of  three  dimensional  dynamics  and  control  of  large  flexible  space 
structures  . 
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5.  CONCLUDING  COMMENTS 


The  analysis  and  design  of  optimal  LQG  digital  shape  and  orientation 
control  for  an  orbiting  shallow  spherical  shell  system  has  been 
studied.  In  designing  for  the  placement  of  the  controller  and  ob* 
server  poles,  it  is  seen  that  the  minimum  modulus  of  the  eigenvalues 
of  the  closed*loop  observer  system  must  be  less  than  the  minimum 
modulus  of  the  eigenvalues  of  the  closed*loop  controller,  so  that 
the  observer  can  provide  for  a  timely  accurate  estimate  of  the  state 
variables  for  the  controller. 

The  problem  of  determining  the  number  and  location  of  the  actuators 
is  also  studied  by  means  of  three  concept  definitions  for  the  degree 
of  controllability.  These  definitions  are  based  on  the  scalar 
measure  of  the  controllability  Grammian  matrix  for  discrete^time 
systems,  and  are  interpreted  both  physically  and  geometrically.  LQG 
transient  responses  for  several  combinations  of  actuator  placements 
verify  the  concept  of  degree  of  controllability  in  selecting  the 
number  and  locations  of  the  actuators. 

This  work  is  also  being  extended  to  the  analysis  and  design  of  opti* 
mal  LQR  digital  shape  and  orientation  control  of  flexible  orbiting 
platform  systems. 

It  is  suggested  that  future  related  research  should  concentrate  on 
the  development  of  practical  design  methods  for  ensuring  the  robust 
control  of  sampled  data  large  flexible  space  systems.  The  main  ob* 
Jective  should  be  to  guarantee  a  certain  minimum  level  of  perform* 
ance  and  stability  robustness  in  the  presence  of  variations  Trom  the 
ideal  design  conditions,  and  under  the  influence  of  unmodelled  (or 
incompletely  modelled)  disturbances.  Specific  applications  to  the 
dynamics  and  control  of  the  orbiting  shallow  spherical  shell  and 
platform  system  could  be  emphasized. 
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APPENDIX 


SEVERAL  PROOFS  FOR  THE  DEGREE  OF  CONTROLLABILITY  PROBLEM 

1.  The  Proof  of  Invariability  of  the  Degree  of  Controllability 
Under  Orthogonal  Linear  Transformation 

It  is  supposed  that  the  linear  discrete^time  system  is  as  follows: 


Xk  =  *xk-i  ♦  ruk-i 

( A  1  ) 

• 

«c 

let  xk  =  Lxk 

( A2 ) 

where 

* 


L  =  orthogonal  linear  transformation 

The  equation  of  the  new  system  after  linear  transformation  is  as 
foil ows : 

xk  =  *xk-M  +  ?  uk-i  (A3) 

where 

*  =  =  lt*l,  l"4  1  r  =  ltt  =  r 

The  controllability  matrix  for  the  system  (A1 )  is 

Pc  =  (r  :  4>r  :  <J»2r  :  ...  4>n-“-1r)  (Ait) 

The  degree  of  controllability  matrix  Wc  is 

wc  =  PcpcT  =  rrT  +  4,ffT*T  +  ...  +  $n-‘1rrT($n-1  )T 
The  controllability  matrix  for  the  system  (A3)  is 

Pc  =  (  r  :  *  r  :  *2  r  :  ...  :  $n,4lr)  =  (L*1  r  :L’‘1  GLL*1  T  :  ...  : 

(L*1  <t>L)  (L*1  4>L)  ...  (  L“  ^  #L  )  L“  ^  T  ) 

=  L‘1(r:$r  :  ...  :  4>n*11D  =  L^Pc 

«c 

4  The  degree  of  controllability  matrix  for  the  system  (A3),  Wc ,  is  as 

follows: 

Wc  =  PCPCT  =  L"1PcpJl  ( A5 ) 

As  we  know,  the  eigenvalues  of  a  matrix  are  invariant  under  a  simi* 
larily  t ransf ormat i on ,  i.e., 

A(WC)  =  A(L'1PcPcTl)  =  A(PcPct)  =  A  ( Wc )  (A6) 
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Therefore,  the  degree  of  controllability  is  also  invariant  under  a 
similarity  transformation,  i.e., 


u i (Wc )  *  u i (Wc )  (i  -  1,  2,  3) 


( A  7 ) 


2.  The  Proof  of  the  General  Properties  (a)  -  (d)  for  the  Degree  of 
Controllability  ui(Wc)  (i  -  1,  2,  3) 

( 1 )  The  case  for  u i 

It  is  evident  that  the  properties  (a)  -  (c)  for  are  true.  What 
we  need  to  prove  is  property  (d).  That  the  property  (d)  is  true  is 
due  to  the  following  facts. 


If  Wc  -  Wc1  + 
U1  (  Wc  )  *= 


wc2»  wc1»  wc2eWc*»  then 

xmin(wc)  -  Min  <WCX,X>  -  Min{  <Wc1X,X> 
|  |  x  |  |  -1  II  x  | |  -  1 

xeRn 

<Wc2X,X>|  >  Min  <Wc1X,X>  +  Min  <Wc2X,X> 
I  |  x  |  |  -  1  I  I  x  |  |  -  1 

xmin(WciX)+Xmin(Wc2)«  Pl(Wc1)+y1(Wc2) 


( A8 ) 


(2)  The  case  for  y2 


Because  the  trace  of  a  positive  definite  symmetrical  matrix  is  invar-* 
iant  under  a  similarity  transformation,  it  13  evident  that  the 
properties  (a)  -  (c)  are  true.  The  proof  of  property  (d)  for  U2  is 
as  follows. 


Based  on  the  definition  of  y2,  what  we  want  to  prove  is  the  follow¬ 
ing  inequality: 


y2(Wc>  =  t?W^T  '  trW^rT  +  trWcT'1  "  U2(Wcl)  +  w2(Wc2)  (A9) 

where  Wc1,  Wc2  e  Wc*  ,  Wc-  Wc1  +  Wc2 
It  is  well  known  that 

Wc"1  -  AdjWc/det  Wc 

then 

n 

trWc-1  -  l  detWc(i)/detWc  ( A 1 0 ) 

i  =  1 

where  detWc(i)  is  the  determinant  of  the  ith  principal  minor  of  the 
matrix  Wc .  Considering  Eq .  (A10),  the  inequality  (A9)  we  want  to 
prove  will  become 
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_ n _  >  n _ _ _  +  _ n 

n  det(Wc1 ( 1 )  +  WC2( i ) )  ~  n  detWci ( i )  n  detWC2( i  ) 

][_,  det(Wc1  *  Wc2)  detwc1  ^  detwc2 


« 


l 


* 


< 


From  the  well  known  Bergstrom  inequality  [R1],  we  have 

det(Wn1+W„?)  >  detW^ i  +  detW„p 

det (Wc1 ( i )*WC2< 1 ) )  detWc i (  i  )  detWc2(i) 

From  ( A 1 2 ) ,  we  have 

det(W„i(i)  ♦  Wn?(l))  <  _ 1 _ 

det (Wc i  +  Wq  2 )  de  t Wn i  +  de  tWn  p 

detwc1 ( 1 )  detWc2( i  ) 

then  we  also  have 

n  det (Wci ( i )+WC2( i ) )  <  n  (  detWc^  detWC2  ^  i 

^  det(Wcl+Wc2)  t  ^  UetWcl  (  i  )  +  detWc2(  i  )  ^ 


i  .e  .  , 


f  n  det(Wci  ( i  )+Wc2(i))^i  >  f  n  (  detwc1 

det  ( Wc  i  +Wc2 )  ]  ~  i./ietWc1(i) 


detWc2 
detWc2( i ) 


j-1)- 


Because  the  detWci /det'rfci  (  i  ) ,  detWC2/detWC2( i  )  are  all  real 
by  applying  the  Minkowski  inequality  [R2] 


(  l 


i-1 


( 


detWc1 
det (Wc1 (  1  ) 


+ 


detW02 
detWj2(  1  ) 


> 


r  y  r _ _ _ 

l1f1  ldet(Wc1(i) 


n  detWc1(i)  n  detWC2(i) 

r  y  - - — i-i  ♦  r  y  - - — i-i 

if,  detwc1  J  if,  detwc2 


Substituting  ( A 1 5 )  into  inequality  ( A 1 4  )  ,  we  have 


r 


n  det  (..C1  (  i  )+Wc2(  ^  )  )  n  detWc1(i) 

y  - , _ 1-1  2  f  y  _  1-1 

1.1  de  t ( Wc i + Wc  2 )  lif1  de  tWc i  1 


n  detwc2( i ) 
detWc2 


( A  1  1  ) 


(  A 1  2 ) 


(  A 1  3  ) 


1  (AH) 

numbers  , 


(  A 1  5  ) 


)-1 

(A16) 
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If  the  two  sides  of  Inequality  (A16}  are  multiplied  by  n,  the  result 
is  just  what  we  need,  (All). 

(3)  The  case  for  113 

It  is  evident  that  the  properties  (a)  ■*  (c)  for  M3  are  true  due 
directly  to  the  definition  of  M3.  Based  on  the  Oppenhelm  inequality 
[ R1 ] ,  we  have 

(de  t  ( Wc  1  +  Wc2))1/n  >  (detWc1)1/n  ♦  (detWc2)1/n  ( A 1  7  ) 

* 

This  is  just  what  we  want  to  show  property  (d)  for  M3,  i.e., 

U3(wc1  +  wc2)  1  ^3(wcl)  +  Wc^wc2^  ( A 1 8 )  » 


t 


A 
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